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Abstract 

Fe-treated biochar and raw biochar produced from macroalgae are effective biosorbents of metalloids and metals, 
respectively. However, the treatment of complex effluents that contain both metalloid and metal contaminants presents a 
challenging scenario. We test a multiple-biosorbent approach to bioremediation using Fe-biochar and biochar to remediate 
both metalloids and metals from the effluent from a coal-fired power station. First, a model was derived from published 
data for this effluent to predict the biosorption of 21 elements by Fe-biochar and biochar. The modelled outputs were then 
used to design biosorption experiments using Fe-biochar and biochar, both simultaneously and in sequence, to treat 
effluent containing multiple contaminants in excess of water quality criteria. The waste water was produced during ash 
disposal at an Australian coal-fired power station. The application of Fe-biochar and biochar, either simultaneously or 
sequentially, resulted in a more comprehensive remediation of metalloids and metals compared to either biosorbent used 
individually. The most effective treatment was the sequential use of Fe-biochar to remove metalloids from the waste water, 
followed by biochar to remove metals. Al, Cd, Cr, Cu, Mn, Ni, Pb, Zn were reduced to the lowest concentration following the 
sequential application of the two biosorbents, and their final concentrations were predicted by the model. Overall, 1 7 of the 
21 elements measured were remediated to, or below, the concentrations that were predicted by the model. Both metalloids 
and metals can be remediated from complex effluent using biosorbents with different characteristics but derived from a 
single feedstock. Furthermore, the extent of remediation can be predicted for similar effluents using additive models. 
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Introduction 

Coal-fired energy generation produces metal-contaminated 
waste water with high concentrations of regulated elements such 
as Al, As, Cd, Mo, Se, V and Zn. The effluent is created when 
water is used to wash residual ash from combustion chambers and 
flue stacks [1]. As these effluents are typically too contaminated to 
be discharged, and current treatment options are restricted by cost 
[2], the effluent from coal-fired power stations is often stored in 
large Ash Dams (ADs). However, high bioavailability and rapid 
bioaccumulation of toxic elements from Ash Dam Water (ADW) 
has resulted in cases of ecotoxicity to vertebrates in ADW 
catchments [3,4] and, therefore, sustainable water treatment 
technologies are required. Biosorption, the use of dried biological 
material to passively bind contaminants from waste water, is an 
option for the treatment of industrial effluents [5,6]. We have 
recently demonstrated that dried freshwater macroalgae are an 
effective feedstock for the production of biosorbents, and native 
species from industrial facilities can be intensively cultivated to 



provide biomass for biosorption [6,7]. This circumvents one of the 
critical constraints to algal-based biosorption, which is sourcing a 
sustainable feedstock for production of biosorbents that does not 
compete with estabhshed markets for cultivated macroalgae [8] . 

Biosorption exploits charge-based interactions between dis- 
solved elements and biosorbents. Negatively charged functional 
groups on the surface of dried macroalgae makes it an effective 
biosorbent with which to remove positively charged metal cations 
from solution and various macroalgal biosorbents have proven 
effective against a range of dissolved metals [9-11]. However, 
existing biosorption research is largely limited to synthetic effluents 
with only one or a few elements targeted for remediation 
[5,7,12,13]. In reality, industrial effluents contain a myriad of 
coexisting elements and the complexities of treating these effluents 
are not replicated in experiments with synthetic solutions [14]. A 
continued focus on the kinetics of metal uptake from synthetic 
effluents by biosorbents, at the expense of empirical data on the 
performance of biosorbents in real-world effluents, has arguably 
limited the application of algal-based biosorption at scale [6, 1 5] . 
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The treatment of complex effluent is difficult because the 
dissolved ions in the effluent have a variety of properties, and 
hence, affinities for biosorbents [6,16,17]. Macroalgae are in 
general effective against a range of dissolved cations and the 
afiinity of dried macroalgae for metals (e.g. Ni^^) can be enhanced 
by converting biomass to biochar through slow pyrolysis [14,18]. 
However, both dried biomass and biochar have low affinities for 
metalloids (e.g. Se04^ ) that co-occur with metals in industrial 
waste water as oxyanions [19-21]. An additional conversion 
procedure, the treatment of biochar with an iron (Fe'^""") solution, 
can enhance the affinit)' of biochar for metalloids [22]. The 
resulting Fe-biochar binds a range of metalloids including Se, As 
and Mo at a rate that is independent of pH but also releases a suite 
of metals (cations), that are natural components of the biomass, 
back into solution in the process [14,22]. While the release of 
metals during deployment of Fe-biochar for the purposes of 
metalloid biosorption may appear counter-productive to the aims 
of treating waste water, sequential application of Fe-biochar to 
remove metalloids and then biochar to remove metals can, in 
principle, treat a complex effluent far more comprehensively than 
if any one biosorbent was applied individually. The application of 
multiple biosorbents to hohstically treat complex waste water 
through biosorption has never been empirically tested. 

Here we test the applicability of a multiple biosorbent treatment 
technique, using Fe-biochar and biochar produced from freshwa- 
ter macroalgae, to treat a complex effluent containing both 
metalloids and metals. We first develop a predictive model for the 
removal of metalloids and metals by Fe-biochar and biochar, 
respectively, from ADW collected from Tarong power station 
(Queensland, Australia). We then apply the model to predict the 
quantity of each biosorbent that is required to maximise the 
removal of metalloids and metals from ADW, accommodating for 
any leaching of contaminants from the biomass during the use of 
Fe-biochar. Finally, we produce Fe-biochar and biochar from 
cultivated green macroalgae Oedogonium, that was originally 
isolated from the AD of Tarong power station, and quantitatively 
compare the efficacy of Fe-biochar and biochar deployed 
simultaneously, and in sequence, to remediate metalloids and 
metals from the ADW. 

Materials and Methods 

Industrial effluent 

This study focused on ADW from Tarong coal-fired power 
station in south-east Queensland, Australia (26.76°S, 151.92°E). 
Tarong is one of Queensland's largest coal-fired power stations 
with a generation capacit)' of 1400 MW and a 46,000 ML AD to 
store waste water from ash disposal processes on site. The ADW 
was sourced directly from the AD and transported to James Cook 
University (JCU), TownsviUe in clean plastic 1000 L Intermediate 
Bulk Containers (IBCs) and stored at ambient temperature in 
12,000 L storage tanks until use. The ADW was collected and 
shipped with the permission and assistance of Stanwell Energy 
Corporation. 

Algal biosorbents production and preparation 

Oedogonium sp. (Genbank: KF606974 [23]) (hereafter Oedogo- 
nium), was used as the feedstock for the biosorbents. Oedogonium 
is a filamentous, freshwater macroalga that is native to Tarong 
AD. The biomass for this study was cultivated in f/2 media in 
2500 L tanks during the austral summer months (January - 
March) in the aquaculture facility on the TownsviUe campus of 
JCU (19.33°S, 146.76°E). Prior to experiments, 2 kg of algae was 
harvested from the tanks and oven dried at 60°C for 48 hours (h). 



The biomass was then converted into biochar by slow pyrolysis 
under conditions previously described [24]. Briefly, Oedogonium 
was suspended in a muffle furnace and purged with N2 gas at 
4.0 L min~' while being heated to 450°C for 1 h. A sample of the 
biochar was converted to Fe-biochar by soaking it in a 5% Fe 
solution (diluted Sigma Aldrich 45% w/v FeClj stock solution) at a 
density of 25 g L"' for 24 h on a shaker plate (100 rpm) at 20°C. 
The Fe-biochar was filtered from the FeCls solution and rinsed 
three times with deionized (DI) water at a rate of 20 ml g~' and 
then dried at 60°C for 48 h. 

Derivation of predictive sorption model 

A model was developed using data from a previous study with 
Tarong ADW [14] (Table SI) to predict the change in 
concentration of 21 elements after the deployment of each 
biosorbent. Most metal and metalloid sorption occurs within the 
first hour of exposure and the most effective sorption occurs when 
tiie pH of die ADW is un-manipulated (pH -7.1) [14]. The 
biosorption data collected under these conditions was used to 
construct the predictive model. Following 1 h of exposure of Fe- 
biochar or biochar to ADW the q-value (the mass of an ion [|j,g] 
adsorbed from, or released into, solution per unit of biosorbent 
[g]), was calculated according to Volesky (2007) [7]. The q-value 
was determined for each combination of biosorbent (Fe-biochar 
and biochar) and element (Al, As, Ba, B, etc.) to derive a model 
that predicted the additive effect of multiple biosorbents, applied 
sequentially or simultaneously, on the remediation of dissolved 
elements within a single effluent: 

Mi= {CxV) - J2[^,xS,j] 

Where, M; is the mass of an element in solution (|ig) following 
treatment number «; C is the initial concentration of the element in 
solution ()j,g L '); V is the volume of effluent solution (L); q,; is the 
q-value of the element for biosorbent, x (|Xg g '); Sx,i is the mass (g) 
of biosorbent, x, added for treatment, i. 

Biosorption experiments 

The aim of the biosorption experiment was to maximise the 
removal of oxyanionic elements (As, Mo and Se) with Fe-biochar, 
followed by the removal of cationic contaminants with secondary, 
or simultaneous, deployment of biochar. As described above, this 
sequential treatment strategy is logical because the targeted 
removal of metalloids by Fe-biochar also contributes some metals 
into solution. This requires a second phase of biosorption with 
biochar to remove existing and leached metals in ADW. 

In total there were 6 biosorption treatments: (1) Fe-biochar 
("FeBC"); (2) biochar ("BC"); (3) sequential Fe-biochar ("FeBC 
FeBC"); (4) sequential biochar ("BC BC"); (5) sequential Fe- 
biochar and biochar ('TeBC BC"); and (6) simultaneous Fe- 
biochar and biochar ("FeBC + BC") (Figure 1). As the aim of the 
biosorption experiment was to achieve remediation of metalloids 
and metals, the primary treatments of interest are the sequential 
and simultaneous application of Fe-biochar and biochar (treat- 
ments 5 and 6 respectively. Figure 1). Treatments 1-4 were 
included as controls to assess the effects of Fe-biochar and biochar 
in isolation for comparison with the sequential and simultaneous 
treatments. Each treatment strategy was repeated 3 times to 
evaluate the model predictions. 

The predictive model was used to select the loading densities of 
each biosorbent for the two main treatments and controls. 
Metalloids were targeted for removal by using a density of Fe- 



PLOS ONE I www.plosone.org 



2 



July 2014 I Volume 9 | Issue 7 | e101309 



Sequential Application of Algal Biosorbents for Bioremediation 



Treatment 1 
(Fe-BC) 



Treatment 2 
(BC) 



Treatment 3 
(Fe-BC Fe-BC) 



Treatment 4 
(BC BC) 



Treatment 5 
(Fe-BC BC) 



Treatment 6 
(Fe-BC + BC) 



4 



ADW 



4 



ADW 



4 



ADW 



4 



ADW 



4 



ADW 



4 



ADW 



Fe-biochar 
(13.7 g/L) 



Biochar 
(60 g/L) 



Fe-biochar (13.7 g/L) 
Biochar (46.3 g/L) 



4 



Anions J., Cations f 



Anions -, Cations [ 



Fe-biochar 




Fe-biochar 




(13.7 g/L) 




(13.7 g/L) 





6 



Anions J.J., Cations fT 



Biochar 




Biochar 




(60 g/L) 




(60 g/L) 





Fe-biochar 




Biochar 




(13.7 g/L) 




(60 g/L) 





Anions -, Cations H 



Anions J,, Cations J, 



Anions ?, Cations ? 



Figure 1. Experimental treatments and predicted changes in anion and cation concentrations in ADW. Arrows signify the direction and 
magnitude of predicted change (up = increase, down = decrease, dash = no change, ?= unknown response). Six biosorption treatments were 
tested: (1) Fe-biochar ("FeBC"); (2) biochar ("BC"); (3) sequential Fe-biochar ("FeBC FeBC"); (4) sequential biochar "BC BC"; (5) sequential Fe- 
biochar and biochar ("FeBC BC" and (6) simultaneous Fe-biochar and biochar ("Fe-BC + BC"). 
doi:1 0.1 371 /journal.pone.01 01 309.g001 



biochar that was predicted by the model to result in removal of 
Mo, the most abundant metalloid, to below ANZECC water 
quality criteria. The result was that Fe-biochar would be used at a 
density of 13.7 g L~' (see further justification in results section 
"Development of treatment scenarios using the predictive model"). A 
limitation of the experiment was that the biochar density predicted 
by the model to remove all metals (1 10 g L ), is not feasible as the 
maximum possible loading density of biochar or Fe-biochar was 
60 g L~'. Therefore, all biochar treatments were loaded at a 
density of 60 g L and all Fe-biochar treatments were loaded at a 
density of 13.7 g L '. The simultaneous "Fe-BC -I- BC" treatment 
included 13.7 g L ' Fe-biochar and 46.3 g L ' biochar to give 
the maximum 60 g L ' of biosorbent (Figure 1). 

Each replicate consisted of a 250 ml plastic beaker with 50 ml 
of ADW and the appropriate mass of biosorbent. The pH of the 
ADW was unaltered (7.06±0.01) for all treatments. Fe-biochar or 
biochar (depending on treatment) was added to the ADW and 
placed in a shaker cabinet (Eppendorf Innova 44R) at 1 00 rpm in 
20°C for 1 h. After 1 h, the biosorbent was separated from ADW 
by centrifugation (7000 rcf, 5 min), followed by two stages of 
filtration (75 |Xm nylon fdter paper, then 0.45 |im syringe filter). 
The filtered solution was transferred to another 250 ml plastic 
beaker for the next treatment (biochar or Fe-biochar depending on 
the scenario). The process of the second sequential treatment was 
identical to the first, except that the mass of biosorbent was 
adjusted to account for any loss of solution during the separation 
process to ensure the same density was applied. After the second 
treatment the solution was processed as described and the water 
samples retained for analysis. 



It is important to note that testing multiple treatments in the 
biosorption experiments provides an opportunity to validate the 
performance of the model when initial elemental profiles differ 
from those under which the model was produced. For example, 
after the first stage of the "BC BC" treatment, the effluent is 
predicted to have low metal, but unchanged metalloid, concen- 
tration. If the model accurately predicts the performance of the 
second BC treatment this indicates the model is robust to 
variations in initial effluent characteristics. Such a finding would 
support the use of the model across diverse effluents with a wide 
range of elemental profiles. 

Elemental analysis 

The concentrations of metals and light metal ions (Al, Ba, Ca, 
Cd, Co, Cr, Cu, Fe, K, Mn, Mg, Na, Ni, Pb, Sr, V, and Zn) and 
metalloids (As, B, Mo, Se) were measured using a Bruker 820-MS 
Inductively Coupled Plasma Mass Spectrometer (ICP-MS; Al, Ba, 
Cd, Co, Cr, Cu, Fe, Mn, Ni, Pb, Sr, V, Zn) or a Varian Liberty 
series II Inductively Coupled Plasma Optical Emissions Spec- 
trometer (ICP-OES; Ca, K, Mg, Na). An external calibration 
strategy was used for both instruments, where a standard solution 
of 0.45 |J,m fdtered ADW was used as the vector to calculate the 
concentration of elements. CoUisional Reaction Interface (CRI) 
was used for As (CRI gas: H2) and V (CRI gas: He), while "^Se 
isotope was used for Se quantification, to eliminate polyatomic 
interferences for these elements. A 1% HCl solution was spiked 
with 1 ppb As, Se and V and measured three times for quality 
control; recovery between 98.5 and 110% indicated no significant 
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interferences. All analyses were conducted at the Advanced 
Analytical Centre at JCU, TownsvUle. 

Data analysis 

Multivariate patterns in biosorption were visualized using non- 
metric Multi-Dimensional Scaling (nMDS) from a Bray-Curtis 
similarity matrix following forth-root transformation. A Permuta- 
tional Multivariate Analysis of Variance (PERMANOVA) was 
conducted in Primer 6.1.14 between the elements and the factor of 
treatment (Figure 1). Univariate analysis took the form of one-way 
Analysis of Variance (ANOVA), with the final concentration for 
each sequential treatment analysed. Data were examined for 
normality and homogeneity of variance using normal-probability 
plots of raw residuals and predicted-residual scatter plots [25]. 
When necessary the data were log-transformed. The PCA and 
ANOVA tests were conducted in Statistica (Ver. 10, Statsoft Inc.). 
The log-predicted final concentrations of each element from the 
biosorption model were plotted against the log-observed values on 
a scatterplot for validation of the model. A trend Une of y = x was 
plotted to indicate where the points would be expected to lie if the 
observed and perfect values were in perfect agreement. Any point 
with a residual value > ± 1 was highlighted as a deviation from the 
predicted concentration. 

Results 

Characteristics of ADW 

Twelve (Al, As, B, Cd, Cr, Cu, Pb, Mn, Mo, Ni, Se, and Zn) of 
the 21 elements measured in the ADW have trigger levels 
established by the Australian and New Zealand Environmental 
Conservation CouncU (ANZECC, Table 1) [26]. Nine of these 
elements (Al, As, B, Cd, Cu, Mo, Ni, Se and Zn) were present in 
untreated ADW at concentrations in excess of the ANZECC 
trigger values for the protection of aquatic Ufe at the 95% level 
(Table 1). These elements are therefore the focus of the following 
results and discussion. 

Development of treatment scenarios using the predictive 
model 

The model was used to predict the amount of Fe-biochar and 
biochar required to achieve comprehensive remediation of 
metalloids and metals to below ANZECC trigger values. Allowing 
for a 15% deviation in removal capacity, the model predicted 
13.7 g L ' of Fe-bioc:har would be sufficient to reduce the 
concentrations of As, Mo and Se in ADW to below the trigger 
values (Figure 2a, c and e). As the Fe-biochar leaches some metals 
into solution [14,22], it was apphed first in the sequential 
treatments so that the subsequent treatment of biochar could 
adsorb both the existing and leached metals from ADW 
(Figure 2b, d, f). The density of biochar required during the 
second treatment for removal of the metals was therefore 
calculated assuming 13.7 g L~' Fe-biochar had been used in the 
first stage of treatment. A density of 1 10 g L ' biochar was 
predicted for the removal of the leached and existing metals. 
Given physical limitations the maximum stocking density that 
could be achieved was 60 g L Consequently, all biochar 
treatments were applied at this dc-nsity (Figure 1). The simulta- 
neous treatments were tailored to result in the maximum removal 
of both metalloids and metals within the physical constraints of the 
biochar system. Therefore, Fe-biochar was applied at a density of 
13.7 g L ' and biochar was added at 46.3 g L giving the 
maximum stocking density of 60 g L~' in treatment 6 ("FeBC + 
BC", Figure 1). The simultaneous treatment ("FeBC +BC", 
Figure 1) was included despite the need for difierent stocking 



densities to determine if a mixed biosorbent consisting of Fe- 
biochar and biochar has potential as a single-step biosorption 
treatment for complex effluents. As biosorption and metal leaching 
occurs relatively rapidly (15 min), it is possible there are rapid 
interactions between leached and bound elements in the 
simultaneous exposure that are not evident in the sequential 
treatments. 

Evaluation of treatment scenarios with biosorption 
experiments 

Treatment of ADW with Fe-biochar and biochar resulted in a 
significandy different final composition of dissolved elements than 
treatment with only Fe-biochar or only biochar (PERMANOVA: 
"treatment" pseudo-F(3_8) = 79.1, p = 0.001). ADW treated with 
Fe-biochar (treatment 1 "FeBC" and treatment 3 "FeBC 
FeBC", Figure 1) had lower concentrations of metalloids (As, Mo 
and Se) than untreated ADW (Figure 3). ADW treated with 
biochar (treatment 2 "BC" and treatment 4 "BC BC", 
Figure 1) had lower alkali, alkaline and transition metal concen- 
trations (Na, Mg, Sr, Ca, Cr, Pb, Fe, Al, Ni, Cu, Zn and Co) than 
untreated ADW (Figure 3). When ADW was treated with both Fe- 
biochar and biochar (treatment 5 "FeBC —> BC" and treatment 6 
"FeBC -I- BC", Figure 1), it had a unique elemental composition 
compared to the other treatments, achieving significandy lower 
concentrations of both metalloids and metals than untreated ADW 
(Figure 3), with the best outcome being from the sequential 
treatment (treatment 5 "FeBC BC") which attained the lowest 
concentrations of metalloids and metals of all the treatment 
scenarios (Table 2). 

a) Sorption by Fe-biochar alone ("FeBC" and "FeBC 
FeBC"). The metalloids (As, Mo and Se) were removed from 
ADW by Fe-biochar (Figure 4a-c; Table 2). Two applications of 
Fe-biochar in sequence (treatment 2 "FeBC FeBC", Figure 1) 
resulted in greater removal of metalloids than any other treatment, 
but also resulted in the highest concentrations of dissolved metals 
being released (Figure 4, Figure SI). For example. Mo decreased 
from 3913 to 43 |xg L"' in ADW (Figure 4b, Table 2), while Zn 
increased from 36 to 2167 |xg L ' (Figure 4e, Table 2). The first 
deployment of Fe-biochar was more efiicient at removing the 
metalloids than the second sequential deployment, which yielded 
lower removal rates of As, Mo and Se per unit Fe-biochar used 
(Figure 4a-c). Therefore, ADW exposed to "FeBC" and "FeBC 

FeBC" treatments was characterised by lower metsdloid and 
higher metal concentrations than untreated ADW (Table 2). 

b) Sorption by biochar alone ("BC" and "BC 
BC"). Metals (Al, Cd, Cu, Mn, Ni, Pb and Zn) were removed 
from ADW by biochar (Figure 4d-f, Table 2). The final concen- 
trations of Al, Cu, Mn, Ni, Pb, and Zn were lowest following the 
sequential biochar treatment ("BC BC") (Figure 4d-e, 
Table 2). Al, Ni and Zn were all decreased to below their trigger 
levels in ADW during the sequential treatment "BC BC" 
(Figure 4d-e, Table 2). As was also reduced in concentration in 
the "BC — > BC" treatment, however, the final concentration was 
higher than the final concentration in the "FeBC — » FeBC" 
treatment (Figure 4a, Table 2). The change in the concentration 
of K was the opposite to the other metals, with substantial leaching 
occurring, resulting in a total increase in tlu' concentration in 
ADW of over 1,800,000 |j,g L"' in tire "BC BC" treatment 
(Figure SI, Table 2). 

c) Sequential and simultaneous treatments ("FeBC 
BC" and "FeBC + BC"). The sequential treatment of ADW 
using Fe-biochar to remove metalloids followed by biochar to 
target the existing and leached metals (treatment 5 "FeBC —> BC") 
resulted in the most comprehensive treatment of the ADW 



PLOS ONE I www.plosone.org 



4 



July 2014 I Volume 9 | Issue 7 j e101309 



Sequential Application of Algal Biosorbents for Bioremediation 
Table 1. Concentration of elements in Ash Dam Water and associated ANZECC trigger values. 



Elem6nt 


ANZECC Trigger (fig L^) 


Initial Concentration [^g L ^ 


( + SE)] 


Aluminium 


55 


89 


(2.1) 


Arsenic 


13 


34 


(0.2) 


Boron 


370 


3823 


(18) 


Cadmium 


0.2 


2.1 


(0.01) 


Chromium 


1 


0.1 


(0.02) 


Copper 


1.4 


2.1 


(0.03) 


Lead 


3.4 


0.03 


(<0.01) 


Manganese 


1900 


0.9 


(0.2) 


Molybdenum 


34 


3913 


(7.2) 


Nicltel 


1 1 


38 


(0.3) 


Selenium 


11 


97 


(0.3) 


Zinc 


8 


36 


(2.0) 


Barium 




100 


(0.3) 


Calcium 




335,333 


\ I Z,\J 1 Z) 


Cobalt 




0.3 


(<0.01) 


Iron 




5.0 


(<0.01) 


Magnesium 




99,800 


(3348) 


Potassium 




44,867 


(2178) 


Sodium 




446,333 


(11,169) 


Strontium 




4080 


(21) 


Vanadium 




982 


(8.1) 



Bold values exceed the ANZECC trigger value for protection of aquatic life at the 95% level. 
doi:10.1371/journal.pone.0101309.t001 



(Table 2). As described above, metalloids were removed by the 
initial application of Fe-biochar while the metals, Al, Cr, Cu, Ni, 
Pb and Zn all leached ofT the biosorbent into solution (Figure 4). 
However, during the subsequent application of biochar, most of 
these metals were adsorbed to below initial concentrations (e.g. 
Figure 4d-e), with the exception of Cu and Ni (Table 2). Of the 
eight ANZECC metals (Al, Cd, Cr, Cu, Mn, Ni, Pb, Zn), four (Al, 
Cd, Pb, Zn) were reduced in ADW to the lowest concentrations 
measured across all treatments, and an additional three (Cu, Ni, 
Mn) were reduced to the second lowest concentrations in ADW in 
the "FeBC BC" treatment (Table 2, Table S2). 

When Fe-biochar and biochar were used simultaneously 
(treatment 6 "FeBC + BC") there were higher final concentrations 
of As, Mo and Se in ADW compared to "FeBC FeBC" and 
"FeBC BC" treatments (Figure 4a-c). B was reduced to its 
lowest concentration in the "FeBC + BC" treatment, however, this 
only constituted a 4% drop from the initial concentration 
(Table 2, Table S2). Similarly, Cr was reduced to the lowest 
concentration following a treatment with "FeBC + BC" of 
1.2±0.04 |j,g L , however, this concentration is still an order of 
magnitude higher than the initial concentration of Cr in ADW, 
0.1 ±0.02 Hg L"' (Table 2, Table S2). Both Cu and Ni had 
equivalent concentrations between the treatments of "FeBC 
BC" and "FeBC + BC" (Table 2, Table S2). An interesting result 
was the response of V, which reduced in concentration following 
sequential and simultaneous treatments, but increased slightly in 
concentration with individual treatments of biochar and Fe- 
biochar (Figure SI). 



Evaluation of the model 

The final concentrations of elements following the "FeBC" and 
FeBC FeBC" treatments were close to those predicted by the 
model (r = 0.97 and 0.96 respectively. Figure 5a-b). The treatment 
of ADW with Fe-biochar delivered greater than expected 
reductions in Mo, but was less effective than predicted for As (in 
both single and sequential treatments) and Se (in sequential 
treatments) (Figure 5a and b). The model was slightly less accurate 
for predicting the concentrations following "BC" and "BC BC" 
treatments, as demonstrated by the weaker correlations between 
predicted and observed concentrations (r = 0.88 and 0.69 respec- 
tively. Figure 5c, d). Despite the weaker correlation, the concen- 
tration of 14 of the 21 elements was accurately predicted for a 
single application of biochar (Figure 5c). The treatment of ADW 
in the "BC" treatment resulted in higher than predicted 
concentrations of B and V and lower than expected concentrations 
of Mn (Figure 5c). The deviations from model prediction were 
compounded in the "BC BC" treatment as 1 1 elements (As, B, 
Ba, Cr, Fe, Mn, Na, Ni, Se, V & Zn) out of the 21 had residual 
values > ± 1 . Na had the greatest deviation from expected 
concentration (Figure 5d; Table 2). B, V and Mn had similar 
deviations from predicted concentrations in the "BC BC" 
treatment (Figure 5d). 

In the "FeBC BC" treatment 17 of the 21 elements were 
remediated close to, or below, predicted concentrations (Figure 6). 
Twelve elements out of 2 1 (Al, As, B, Co, Cr, Cu, Fe, Mo, Pb, Se, 
V & Zn) had residuals > ± 1 (Figure 6). However, of these most 
were removed from ADW more effectively than was predicted by 
the model. As, B, Se and V had higher concentrations than 
predicted after the "FeBC BC" treatment (Figure 6). Despite 
the high rate of leaching of Fe from the Fe-biochar treatment, the 
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Figure 2. Modelled changes in concentrations of (a) As, (b) Al, (c) Mo, (d) Ni, (e) Se and (f) Zn In ADW treated with Fe-biochar and 
biochar (solid and dashed lines respectively). ANZECC trigger levels and stocking densities which result in reduction below the ANZECC 
concentration indicated by horizontal grey line and dashed area, respectively. The initial concentration of Al, Ni and Zn for the biochar treatment in 
plots (b), (d) and (f) are as predicted following application of 13.7 g Fe-biochar. The grey shaded area indicates biosorbent densities above the 
physical limit for experimentation (60 g L~^). 
doi:1 0.1 371/journal.pone.01 01 309.g002 



subsequent biochar treatment removed all of this Fe, resulting in a 
concentration (5 |ig L ) much lower than was predicted 
(39,451 ng L"') (Figure 6; Table 2). 

The predictive model is included as supporting information to 
this publication, allowing users to insert metal concentrations and 
desired efiluents to see predicted treatment using the sequential 
approach (Table S3; instructions and assumptions available in the 
table caption). The elements are treated independendy in the 
model with the assumption that the q-value for each element is 



independent of initial concentration. The model also assumes the 
effects of Fe-biochar and biochar use are additive. Therefore, the 
model output will show linear increases in the requisite Fe-biochar 
stocking densities as the initial concentration of each target 
metalloid in solution increases. The amount of biochar required in 
the second step of the sequential treatment will also increase 
linearly with the stocking density of Fe-biochar used. Instructions 
for use are supplied along with Table S3. 
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Figure 3. Non-metric Multi-Dimensional Scaling (nMDS) demonstrating the sorption of elements by Fe-biochar, biochar, and a 
combination of the two biosorbents in single and sequential treatments. (A) nMDS (Stress <0.01) with cluster analysis superimposed. Open 
and closed squares and circles represent single and sequential applications of Fe-biochar and biochar, respectively. Triangles represent sequential 
application of Fe-biochar followed by biochar and diamonds represent simultaneous application of Fe-biochar and biochar. (B) nIVlDS with vectors 
superimposed, the direction and length of indicate the strength of correlation with the treatment. 
doi:1 0.1 371/journal.pone.01 01 309.g003 




Figure 4. Change in solution concentration of (a) As, (b) Mo, (c) Se, (d) Al, (e) Zn, and (f) Cd In ADW after sequential exposure to Fe- 
biochar and biochar. The "BC", "FeBC" and "FeBC BC" treatments are represented as dotted, dashed and solid lines, respectively. The "FeBC + 
BC" treatment is represented by the circle in each panel. The "FeBC + BC" result is placed under treatment 2 to compare with the final concentrations 
of the other treatments. ANZECC trigger level is represented by a horizontal grey line. 
doi:1 0.1 371/journal.pone.01 01 309.g004 
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Figure 5. Scatterplots of predicted vs. observed concentrations of elements in ADW (|ig L^) following exposure to (a) "FeBC", (b) 
"FeBC —>■ FeBC", (c) "BC", and (d) "BC —>■ BC" treatments. The data have been log transformed. The plotted line shows the expected 
relationship if all observed values are identical to predicted values. Points falling above the line were overestimated by the model (higher than 
expected final concentrations). Points falling below the line were underestimated by the model {lower than expected final concentrations). Asterisks 
indicate elements with a residual of greater than ±1. 
doi:1 0.1 371/journal.pone.01 01 309.g005 



Discussion 

We have demonstrated that the sequential apphcation of 
multiple biosorbents provides a comprehensive treatment of a 
complex effluent. The combination of targeted remediation of 
metalloids by Fe-biochar and metals by biochar resulted in a 



greater number of elements being treated compared to the 
application of either biosorbent independently. Furthermore, the 
reduction of metalloids followed by metals, b)' tlu' sequential 
application of Fe-biochar then biochar, demonstrates that the 
effects of these biosorbents are additive. Our model predicted the 
concentration of most elements following treatments with Fe- 



PLOS ONE I www.plosone.org 



10 



July 2014 I Volume 9 | Issue 7 | e101309 



Sequential Application of Algal Biosorbents for Bioremediation 



FeBiochar -> Biochar 



16 
14-1 
12 



r = 0.71, p < 0.001 




-2 0 



8 10 12 14 16 



Log Predicted Concentration (ug L^) 

Figure 6. Scatterplot of predicted vs. observed concentrations of elements in ADW (^g L~^) following the "FeBC —> BC" treatment. 

Data have been log transformed. The plotted line shows the expected relationship if all observed values are identical to predicted values. Points 
falling above the line were overestimated by the model (higher than expected final concentrations). Points falling below the line were 
underestimated by the model (lower than expected final concentrations). Asterisks indicate elements with a residual of greater than ±1. 
doi:l 0.1 371/journal.pone.01 01 309.g006 



biochar, although the model was slightly less efiective at predicting 
the concentrations of elements when biochar was involved in the 
sequence. Regardless, the selectivity of the biosorbents was 
consistent and accurately modelled across the six water types 
(i.e. ADW treated by Fe-biochar, biochar, and combinations 
thereof) in which the initial element composition varied greatly. 
This provides the basis to predict the efficacy of these treatments 



(Fe-biochar and biot:har in any order of application) in remedi- 
ating the eflluents from other industries. 

Treatments that combined Fe-biochar and biochar ("FeBC 

BC" and "FeBC + BC") achieved the most comprehensive 
remediation of the greatest number of elements. However, 
s(;quential treatments with Fe-biochar ("FeBC FeBC") and 
biochar ("BC BC") were more effective at removing their target 
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constituents (metalloids and metals, respectively). In the same way 
that the removal of a given element will be greater in isolation than 
when that element is present in a complex mixture, removing 
multiple elements from an effluent occurs at the expense of the 
sorption capacity of the biosorbent for each individual element 
[27]. Although treatment of a wastewater with a sccjuc-ntial 
application of macroalgal biosorbents has not been done 
previously, it has proven successful with other biosorbents. For 
example, in an analogous experiment, a sequential approach was 
taken to the treatment of wastewater using anaerobic fermented 
biological sludge, first using raw sludge to remove metals, then 
using sludge that had been pre-treated with a cation detergent 
(either Al''^ or Fe'^"*} to remove the metalloids [28] . We found no 
benefit in terms of efficacy for simultaneous application of the 
biosorbents, and some negatives regarding the loading ratios per 
unit volume, meaning water should be treated in two-step process 
of Fe-biochar followed by biochar. 

Predicting the concentrations of elements following treatments 
of Fe-biochar and biochar with the model was very successful, with 
all of the sequential and single treatment scenarios resulting in at 
least a moderately strong correlation between predicted and 
observed values. ModeUing has been used extensively in biosorp- 
tion studies, however often only to predict and quantify the 
mechanism by which the element is being absorbed onto the 
biosorbent using Langmuir and Freundlich isotherms 
[12,16,29,30]. Our empirical model performed particularly well 
when considering that the starting concentrations of most elements 
were significantly different to those in the effluent that was used to 
derive the model, and that the water quality after each treatment 
was unique, spanning 6 t^pes of water treated with different 
combinations of biosorbents, resulting in a broad range of initial 
concentrations (Table 2). In this context, the model is as a robust 
tool for predicting the treatment of a broad range of waste water 
profiles by biosorbents derived from Oedogonium, and as a tool to 
develop working treatment strategies for a variety of waste water 
sources with known elemental compositions. 

The model presented here intrinsically has some assumptions 
and limitations. The model is based on q-values from single 
deployment of the biosorbents in untreated ADW. The key 
assumption of the model is that there is a linear relationship 
between the extent of biosorption and biosorbent density that is 
independent of the initial concentration of each element in 
solution [16]. The model was slightly more effective at predicting 
final concentrations of elements after a single deployment of Fe- 
biochar or biochar than after two deployments. This suggests that 
the q-value is not entirely independent of the initial element 
concentrations because the biosorbents tend to r(;m()\'(; sliglitly less 
of each element than predicted when initial concentrations are 
lower than the conditions under which the model was derived 
[31,32]. However, the magnitude of difference between predicted 
and observed final concentrations are small despite the large 
variation in initial conditions across the six treatments. This 
reiterates the importance of assessing the efficacy of biosorbents 
under realistic initial conditions, rather than deriving models of 
biosorbent kinetics in single-element system with inordinately high 
initial metal concentrations. Regardless, our approach to deter- 
mining q-values of biosorbents for multiple elements in a real- 
world effluent provides a robust metric for use in biosorption 
modelling and with additional testing under different scenarios, 
further accuracy and confidence limits may be applied. 

Conclusions 

There is substantial literature on macroalgal biosorption but the 
majority of existing data has assessed elements in isolation, while 
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an industrial effluent typically contains a multitude of co-existing, 
interacting elements. Given the large quantities of biomass that 
win inevitably be required to treat industrial effluents, it will also 
be necessary to select and cultivate biomass locally for the express 
purposes of bioremediation. We have demonstrated that different 
types of biosorbents (Fe-biochar and biochar) can be produced 
from a single feedstock that is native to an industrial facility (the 
freshwater macroalga Oedogonium, see also [18,33]) and used to 
sequentially treat metalloids and metals from a complex effluent. 
Through the use of a predictive model we demonstrate that 
sequential biosorption is largely predictable and consistent across a 
wide range of initial conditions. Targeting metalloids that do not 
have a natural affinity for biochar requires Fe-treatment, which 
inherentiy results in a suite of additional metals being released into 
solution. This metal leaching can in turn be addressed through the 
use of biochar as a final treatment. Therefore, while biosorption 
has been widely cited as a sustainable and cost-effective means of 
treating waste waters, our data clearly show that the reality is far 
more complex than is typically acknowledged, but remains 
achievable. 

Supporting Information 

Figure SI Change in solution concentration of (a) 
potassium, (b) manganese, and (c) vanadium following 
sequential exposure to Fe-biochar and biochar. Sequential 
exposure of the ADW to biochar, Fe-biochar, and Fe-biochar 
followed by biochar and represented as dotted, dashed and solid 
lines, respectively. Exposure of ADW to both Fe-biochar and 
biochar simultaneously is represented as a round dot. Simulta- 
neous exposure only had one application yet is placed under 
treatment 2 to compare with the final concentrations of the other 
treatments. ANZECC trigger level represented by a horizontal 
grey line. Error bars show standard error. 

fDocx;) 

Table SI Change in dissolved elemental concentration 
((ig L~') in ADW following treatment with biochar and 
Fe-biochar for 1 h at a solution pH of 7.1, and the 
derived q-values (ng g Data were collected during 
experiments for KidgeU et al (in press). 
(DOCX) 

Table S2 Summary table for ANOVA tests run on each 
of the 21 elements investigated. One-way analysis of 
variance tests were run on final elemental concentration with the 
factor of Treatment. Type III sum of squares was used. AU tests 
met the assumption of homogeneity of variance, normality of 
residuals and independence. Transformation of the data were 
required for some elements, the transformation applied is listed 
next to the tide. Factors in bold indicate significance under alpha 
of 0.05. 
(DOCX) 

Table S3 Predictive biosorption model for treatment of 
industrial waste water through sequential Fe-biochar 
and Biochar applications. Instructions to users. 1. Cells 
that can be edited by the user are shaded orange. 2. Input initial 

dissolved concentrations of elements and the volume of effluent 
requiring treatment in the orange shaded cells at Step 1. The 
model will calculate the recommended Fe-biochar stocking density 
which is displayed at Step 2. 3. Select the maximum Fe-biochar 
stocking density for elements of interest from Step 2 and enter it in 
the orange shaded cell at Step 3. The model will calculate the 
stocking density of biochar required to remove leached and 
existing metals in the effluent after Fe-biochar deployment which 



PLOS ONE I www.plosone.org 



12 



July 2014 I Volume 9 | Issue 7 | e101309 



Sequential Application of Algal Biosorbents for Bioremediation 



is played at Step 4. 4. Select the maximum biochar stocking 
density for elements of interest from Step 4 and enter it in the 
orange shaded cell at Step 5. The model will calculate the 
predicted final concentrations of metalloids and metals after the 
sequential application of Fe-biochar and biochar in the efiluent at 
Step 6. 

pa.sx) 
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